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I Abstract. A summaryispresentedofthemostrecentmatrixelementsformassless 2 ^2 scattering 

processes calculated at two loops in QCD perturbation theory together with a brief review on the 
calculational methods and techniques used. 



INTRODUCTION 



The description of high energy processes at hadron colliders together with the study 
of the environment these events are emerging into, will rely heavily on the predictive 
■ power of the theory of Quantum Chromo-Dynamics (QCD). It is well known that QCD 

studies enable a better understanding of the color flow dynamics, which in turn facilitates 



^ I improved predictions on observables such as pj, distributions and production rates. 

Also, QCD studies have a direct impact on theoretical predictions for signals and their 
Oh. backgrounds, which are a key element in the quest for precision physics at the highest 

r-| ! energies. A plethora of data on multi-particle final states will soon become available 

from high-energy collider runs and the comparisons of jet observables with theoretical 
- predictions will be vital. In fact, missing higher order QCD theoretical predictions 

rS I for such observables may be large and important. Therefore a systematic approach to 

[ perform these complex calculations, which usually involve the analysis of multi-leg 

and/or multi-loop amplitudes, becomes imperative. In recent years, new techniques and 
integral manipulation methods have been designed to achieve outstanding analytical 
results in this area. It is the purpose of this talk to review briefly the latest results on 
two-loop QCD corrections for massless 2 — > 2 parton scattering processes and some of 
the techniques used. This is not an exhaustive review, so more details can be found in 
the references provided and the ones therein. 

There are many important reasons why one might consider a next-to-next to leading 
order (NNLO) calculation[|T]] and to mention a few. 

Reduced scale dependence The theoretical prediction for any observable T should be 
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independent of the renormalisation scale /i. But, if we perform a fixed order calcu- 
lation a scale dependence is introduced. The change due to the variation of renor- 
malisation scale is formally one order higher than the one at which the theoretical 
prediction is given. Often, when giving a NLO result, a variation of jl around some 
hard scale is used as a tool to give the uncertainty in uncalculated higher orders but, 
this is an estimate of their size. The theoretical prediction can be improved with a 
complete NNLO result. 
Improved jet description If we consider higher order corrections we automatically 
improve the matching between theoretically and experimentally defined jets. At 
leading order we are modeling jets with single parton emission, but at higher orders 
the phase space available is extended so that at NNLO up to three partons can 
combine to form a jet. Corrections beyond LO involve a better description of soft 
gluon radiation within the jet and this may provide a more accurate picture of its 
shape and structure. 

Reduced power corrections When comparisons of NLO predictions with experimental 
data are made, the need for power corrections (terms of ^(1/2")) arises. The 
structure of these corrections can be motivated theoretically but they are always 
fitted to experiment. One might expect that higher order corrections may play a 
role in reducing the size of these power corrections. 



Partonic cross sections beyond NLO 

The description of hard scattering processes at future hadron colliders requires the 
study of the factorized structure of a cross section for processes with quarks and gluons 
in the initial state. Up to power corrections, the inclusive factorized cross section can be 
written as 

(7(P1,P2) r^Y^j dxidx2fi/i{xi,^j)fj^2(^i,nj)dij{p^,p2,as{^^),s/^^,s/^j) (1) 



where = x^Pf is the momenta of the partons that initiate the hard scattering and /^^^ 

are the parton distribution functions (pdf 's). The non-perturbative effects are comprised 
in these distribution functions (and into fragmentation functions, in the exclusive case), 
a is the hard scattering partonic cross section which describes the interaction of partons 
{ij} that arose from hadrons {1,2} and can be calculated perturbatively, so that at 
NNLO for 2-particle production it can be written as 



^2 jet- 
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(2) 



where [ ]„ indicates the number of particles in the final state with d^n the corresponding 
phase space and ^^'^ the i-th order scattering amplitude. After renormalisation, each 



TABLE 1. Number of Feynman diagrams for different 2^2 scattering processes 



Number of Diagrams 



Process 


Tree 


One loop 


Two loops 


Process 


Tree 


One loop 


Two loops 


gg^gg 


4 


-80 


- 1700 


qiq\ ^qih 


1 


- 10 


-200 


qq^gg 


3 


-30 


-600 


qq^gy 


2 


- 10 


-300 


qq^qq 


2 


-20 


-200 


qq^YY 


2 


- 10 


- 140 



+ crossed processes 



of the integrals in Eq. (||) is ultra-violet (UV) finite but infra-red (IR) divergent which 
manifests itself as poles in e (we adopt D = A — 2e as a dimensional regulator)^]. 
The integration over phase space and the cancellation of poles requires the study of 
the kinematical regions where the additional radiated particles become unresolved^. 
Finally, pdf's and their evolution are needed at an accuracy that matches that of the 
matrix element calculation and for NNLO some great developments have already taken 
place®. 



TWO-LOOP INTEGRALS 

The number of diagrams involved in the complete two-loop and one-loop matrix el- 
ements needed for the NNLO contribution to inclusive jet cross sections and photo- 
production at hadron colliders, are shown in Table |I] [|^]. The types of integrals arising 
in these matrix-element calculations include integrals with scalar numerators 

where fc- = k,i are the loop momenta, p ■ are the external momenta and ^ is a scalar 
function. Here, the massless propagators are denoted by A- ~ (^, ± p+i£)^ where p 
can be any of the loop or external momenta. There are other terms that would typically 
look like the one shown in Eq. (|]), except now the function ^ ^^^"'{k,t} is a ten- 
sor that can depend on the loop-momenta of the system and/or other tensors (such as 
the metric tensor g^"^). Furthermore, the denominator of an integral provides the mo- 
mentum flow in the graph and gives a complete description of its skeleton or topology 
which can be planar or non-planar. Within these two categories we will identify sub- 
topologies or pinchings, depending on whether or not a particular subset of propagators 
is absent. As technical aids, we use a couple of auxiliary general diagrams shown in 
Fig. |l], where each of the propagators is labeled by an integer and carries a specific mo- 
mentum that fulfills conservation of momenta throughout. This auxiliary representation 
allows us to manipulate hundreds of integrals with the minimum amount of informa- 
tion and without compromising the accuracy of the integral description. So to refer to a 



The singularities are guaranteed to cancel for sufficiently inclusive physical quantities[^. 
See for example, analytical!^ and numerical[Q] appUcations. 



Ai = e 

A2 = {k+pif 

A3 = {k+pi +p2f 

A4 = {k+Pi+P2 +P3f 


A5 = P 

A, = {£ + pif 

At = ie + pi +P2f 

As =(£ + pi+p2 +P3f 


Aa = ik- If 


Aw = {k - £- Pi - P2 - Psf 



ft P3 P3 

FIGURE 1. Auxiliary diagrams and table with description of propagators 



Identification 


Diagram 


Identification 


Diagi 


am 


SUNC = /^[0, 1,0, 0,0, 0,0, 1,1,0] 
TRI = /"[1, 0,1, 0,0, 0,0, 1,1,0] 
GLASS = 0, 1, 0, 1, 0, 1, 0, 0, 0] 
XTRI = J^[l, 0, 1, 0, 0, 0,1,1, 1, 1] 
CBOX = I^[l, 1, 0, 0, 0, 0,1,1, 1, 0] 




ABOX = I°[l, 1,1,0, 0, 0, 0, 1,1,0] 
PBOXl = J"[l, 1,1, 0,1, 0,1, 1,1,0] 

PB0X2 = I°[l, 1, 1, -1, 1,0, 1, 1, 1, 0] 
XBOXl = J^[l, 1, 1, 0, 0, 0,1,1, 1, 1] 

XB0X2 = J"[1,1,1,0, -1,0, 1,1,1,1] 
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FIGURE 2. IVIaster integrals needed for the two-loop massless 2 —^ 2 parton hard scattering matrix 
evaluation 



general D-dimensional integral with 10 propagators raised to arbitrary powers, we use 
/^[Vj, V2, V3, V4, V5, Vg, V7, Vg, Vg, Vjq]. Any planar or non-planar topology is represented 
by eliminating, increasing or decreasing the values of v^. It turns out that the two-loop 
QCD matrix-element calculation for 2^2 scattering processes with massless external 
states, requires the explicit analytic evaluation of the 10 master integrals (MI) shown in 
Fig. 1^. Any other integral that contains more powers in the propagators or higher dimen- 
sions can be written as a linear combination of the MI and below we will explore briefly 
how this can be done. First though, we will review the methods that are used to evaluate 
these MI explicitly. 



Explicit loop integration 

The strategically important analytic expressions for all MI have been provided in 
Refs.[0, ^ ^ as series in e = (4 — D)/2. Usually, Feynman or Schwinger parameter- 
izations (see for example [[T(]|]) are used, whereby the propagators of the loop integrand 
are expressed in terms of integrations for real parameters over a particular range. In 
both of these approaches, the loop momenta integration can be done easily and the re- 
maining integrations over the parameters are doable for reasonably sized topologies. 
Sometimes, using these parametric representations does not leave an integral that can 
be solved easily, even more so when the number of loops, external legs and kinematical 
scales increase. Different techniques arose due to this situation, among others 



• The Mellin-Bames method is based on the representation for a sum to some power. 



as a contour integral over a complex variable and the integration is then performed 
on straight contour lines parallel to the imaginary axis. After closing a contour, the 
result is the sum of all enclosed residues that may be expressed as a hypergeometric 
series. This has been used to obtain vital results for the two-loop boxes[^ [Tl]]. 

• The Negative Dimensions technique consists of rewriting the integral over the pa- 
rameters by introducing new ones through a multinomial expansion. Many con- 
ditions have to be satisfied among the parameters, which leads to the restriction: 
D must be a negative integer. Some results for MI have been obtained using this 
method[|ll. 

• Numerical strategies are used, where any loop integral is stripped analytically of 
its IR singularities so that the finite integrals can be performed numerically. These 
methods can also be used to tackle the phase space integrations in the calculation 
of observableslQ]. 

• The analytic evaluation of MI can also be carried out without explicit integration 
over loop momenta by deriving differential equations for MI in internal propagator 
masses or in external momenta. The equations can be solved with appropriate 
boundary conditions. This approach has been widely used to evaluate MI, as was 
done recently for e+e^ 3 



Loop integral reduction 

This approach produces an environment in which complex topologies with high 
powers on the propagators can be reduced down to integrals that can be solved with the 
methods described in the previous section. This reduction can be achieved using systems 
of equations that stem from Integration by Parts identities (IBP)[|T^ and exploiting the 
Lorentz invariance[^ of Feynman integrals. Let us review briefly how these identities 
are generated and help reduce the complexity of an integral. Consider a general two-loop 
integral 
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The idea behind IBP is to generate relations between loop integrals through a total 
derivative with vanishing surface terms, expressed as the following identity 



d'^k 

~ w 
in 2 



° dk^ 



in 2 



0, 



(5) 



where fc- = fc, £ are the loop-momenta and Y can be any internal or external momenta 
involved in the loop integration. Executing the derivative on all possible choices for 
will generate a set of relations^ between integrals with dot products in the numerator. 



For a graph with m loops and n independent external momenta, we can generate m{m + n) identities. 



These dot products can be rewritten in terms of linear combinations of propagators, by 
means of relations such as 

2{k + g).{k + h) = {k + gf + {k + hf-{g-h)\ (6) 

With this simple step we can rewrite all the contents of the numerator in terms of 
propagators that may or may not be part of the denominator. We then say that the 
numerator is reducible if we can cancel it through and irreducible otherwise. In most 
cases, we can exchange the problem of calculating the original integral, for the problem 
of calculating a set of simpler integrals. We can imagine studying all the IBP identities 
for the topologies involved in a particular matrix-element calculation and applying the 
reduction procedure iteratively. Then we can assemble an algorithm that takes any 
integral and expresses it in terms of a few MI. This is precisely what we will review 
next. 



MATRIX ELEMENT EVALUATION AND RESULTS 

The ideas exposed above and extensions to them have been used in the past few years 
to tackle two-loop calculations. The way these calculations are put together varies, 
for example the chart in Fig. ^ shows some basic steps we followed to evaluate the 
two-loop QCD matrix-elements for 2 ^ 2 scattering. Steps 1 through 3 are mainly 
related to the automatic generation of diagrams using QGRAF[p31| and the tensor algebra 
manipulation using FQRM[[T^]. By step 4, we have the matrix element written as a linear 
combination of many different scalar and tensor two-loop integrals. These integrals can 
be rewritten in terms of MI (step 6), using a reduction system that was generated using 
MAPLE (step 5). The final step is the input of e-expansions for each MI, so that we get 



n=0 ^" 

where {s,t, u} are the usual Mandelstam variables. The functions /„ depend on the color 
factors {Np^Cp^Cj^.Tj^ and logarithms of ratios of the kinematic variables. The finite 
piece /q can contain up to 6 polylogarithms (Li„ with n = {2,3,4}) of ratios in the 
kinematical variables, and the Riemann zeta function ^(n) (with n = {2, 3}). 

Several checks have served to verify the results obtained with this approach such 
as comparisons with QED processes [[T5|] and also the study of the analytic structure 
of QCD amplitudes in the limit of forward and backward scattering together with the 
confirmation of the two-loop gluon Regge trajectoryJT^]. The fact that the singular 
analytical structure in Eq. ([7|)can be obtained independently, is also a powerful check 
on the calculation[|^ . CataniQlTl] proposed the structure of the 1/e" V n = {2, 3, 4} poles 
together with the color uncorrelated structure of the 1/e pole. For a while, the origins 
of this proposal did not exist. However, in the work with Sterman[[T8|], we describe how 
the factorization properties of loop amplitudes lead to the exponentiation of double and 
single poles at each order in perturbation theory. The poles can then be assembled 
in terms of universal functions associated with the external partons. This formalism 
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QGRAF 
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in/out states i 
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(STEP 3)-»- 


FORM 


Feynman rules] 
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MAPLE 
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IBP/LI identities 
Solve system 
of equations 

J 







List of all integrals 

in terms of MI : 

int(l,2,l,0,l,0,l,l,l,s,t) = 
a HE +b -e- 
c TTT + d -d + ■ 



-«- (STEP 5) 



(step 4) -»- 



Sum of all tree x loop diagrams 
in terms of integrals : 
(.■Vl(0)|.'V((2)) = 

a int(l,2,l,0,l,0,l,l,l,s,t) 
+ b int(l,0, 1,0, 1,-1, l,l,l,t,u) 



e expansions 





Sum of all tree x loop diagrams 




in terms of MI : 










a HE +b -e- 




c TTT + d -dr + ••• 



-«- (STEP 6) 



fi ~ fi(ln,Li2,Li3,Lii) 



^STEP 7) 



FIGURE 3. Basic steps to follow in the evaluation of matrix elements 



provides a way to generate the complete pole structure for multi-loop amplitudes at 
two loops and beyond. 



CONCLUSION AND OUTLOOK 

As reviewed, the past few years have seen a breakthrough in multi-loop integration 
technology and outstanding progress in the calculation of two-loop matrix elements in 
QCD. Much work remains to be done to have NNLO Monte Carlo numerical estimates 
but, looking back at the results we now have, it seems they will be available soon for the 
first basic scattering processes. This will enable an improved description of high-energy 
QCD phenomena. 
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